On the origin of the hump structure in the in-plane optical conductivity of high Tc 

cuprates based on a SU(2) slave-boson theory 
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An improved version of SU(2) slave-boson approach is applied to study the in-plane optical conduc- 
tivity of the two dimensional systems of high Tc cuprates. We investigate the role of fluctuations 
of both the phase and amplitude of order parameters on the (Drude) peak-dip-hump structure 
in the in-plane conductivity as a function of hole doping concentration and temperature. The 
mid-infrared(MIR) hump in the in-plane optical conductivity is shown to originate from the antifer- 
romagnetic spin fluctuations of short range(the amplitude fluctuations of spin singlet pairing order 
parameters), which is consistent with our previous U(l) study. However the inclusion of both the 
phase and amplitude fluctuations is shown to substantially improve the qualitative feature of the 
optical conductivity by showing substantially reduced Drude peak widths for entire doping range. 
Both the shift of the hump position to lower frequency and the growth of the hump peak height 
with increasing hole concentration is shown to be consistent with observations. 



PACS numbers: 74.20.Mn, 74.25.Fy, 74.25. Gz, 74.25.-q 



I. INTRODUCTION 

High Tc superconductors are believed to be the sys- 
tems of strongly correlated electrons essentially in two 
space dimension. This strong correlation results in the 
peak-dip-hump structure of the optical conductivity a{Lo) 
and linear frequency dependence of the scattering rate 
l/r(u;) [1-6]. Both the appearance of the hump(which 
occurs at w « IOOOcto-^ in YBCO [3,4] and Bi2212 
[6]) and the linear frequency dependence of the scat- 
tering rate(l/r(ci;) ~ w) indicate strong deviations from 
the Drude model prediction of conventional Fermi-liquid. 
Various theories have been proposed to explain these 
non-Fermi-liquid(NFL) like behavior in the charge dy- 
namics (optical conductivity) of the high Tc cuprates. 
Using the nearly antiferromagnetic Fermi-liquid theory, 
Stojkovic and Pines [7] reported a study of normal state 
optical conductivity for the optimally doped and over- 
doped regions, but not for the underdoped regions. They 
showed that the highly anisotropic scattering rate in dif- 
ferent regions of the Brillouin zone leads to an aver- 
age relaxation rate of the marginal Fermi-liquid form by 
showing 1/t(ljj) ^ uj. Their computed optical conductiv- 
ity agreed well with experimental data for the normal 
state of an optimally doped sample. Using the spin- 
fermion model [8,9] and spin susceptibility parameters 
obtained from inelastic neutron scattering(INS) and nu- 
clear magnetic resonance(NMR) measurements, Munzar 
et al. calculated the in-plane optical conductivity of op- 
timally doped YBCO [10]. They examined the peak- 
dip-hump structure only at optimal doping and showed 
a good agreement with observation. From the computed 
self energy they showed that the hump is originated from 
the hot spot and the Drude peak from the cold spot in 



the Brillouin zone. Haslinger et al. reported optical 
conductivities cr(w) of optimally doped cuprates in the 
normal state by allowing coupling between the fermions 
and the bosonic spin fluctuations [11]. They found that 
the width of the peak in spectral function Ak(i^) scales 
linearly with uj in both hot and cold spots and a{u!) is 
inversely linear in uj up to very high frequencies. Be- 
sides these spin fluctuation theories mentioned here, by 
using a phenomenological form of the charge collective 
mode(CM) propagator and empirical parameter values 
Caprara et al. showed that charge-ordering instability 
is responsible for the hump in the optical conductivity 
of high Tc cuprates [12] limited only to the overdoped 
region. The various above theories are reported in the 
limited range of hole doping resorting to empirical pa- 
rameters by fitting INS or angle resolved photoemission 
spectroscopy(ARPES) measurements. Thus the applica- 
bility of these theories may be limited only to a certain 
range of hole doping and temperature. It is, thus, highly 
desirable to use a theory which does not depend on any 
empirical parameters for the entire range of hole dop- 
ing and temperature encompassing both the pseudogap 
phase and the superconducting phase. 

Earher Lee and Salk [13] proposed a U(l) and SU(2) 
slave-boson theories and showed an arch-shaped Tc line 
in agreement with observations. Their theory is different 
from other previous slave-boson theories [14-17] in that 
the Heisenberg term in the t-J Hamiltonian contains the 
contribution of coupling between the spin and charge de- 
grees of freedom(the spinon pair and holon pair orders). 
Recently, using the same U(l) slave-boson theory, we [18] 
were able to explain the peak-dip-hump structure of the 
observed optical conductivity, by showing that the hump 
is caused by the presence of spinon pairing order formed 
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from the hot spot in the Brillouin zone. In this paper, by 
using the SU(2) theory we report a rigorous examination 
on the origin of peak-dip-hump structures over the cnitire 
range of temperature(0 < T < Tc, Tc < T < T* and 
T > T*) and the entire(underdoped, optimally doped 
and overdoped) range of hole doping. In addition we 
discuss how the inclusion of the order parameter phase 
fluctuations markedly improve the predicted in-plane op- 
tical conductivity over the previous results [18]. 



II. THEORY 

A. U(l) AND SU(2) SLAVE-BOSON THEORIES 
OF t-J HAMILTONIAN 

In order to bring forth the essential physical points 

we briefly review the U(l) and SU(2) slave-boson the- 
ory and point out the importance of coupling between 
the charge and spin degrees of freedom in determining 
the stiTicture of the optical conductivity [13]. The U(l) 
formulation is readily done by rewriting the electron op- 
erator as a composite of spinon(/) and holon(6) opera- 
tors, Cicr = fiabl with the single occupancy constraint, 

J2cr fiafi'7~^b\bi = Then we obtain the partition func- 
tion, 



= j VfVbVXe~-fo 



(1) 



with C = E^{J:afl9rf^a+bl^rb?) + Ht-J + 

^liY^a ficrfto- + blbi - 1) where Xi is the Lagrange 
multiplier field to enforce the single occupancy constraint 
and Ht-j, the U(l) slave-boson representation of the t-J 
Hamiltonian, 

Ht-J = -t (flfj^blbi + c.c.) 

<i,j>,a 



<i,j> 



(2) 



It is noted that the above Heisenberg interaction 
term(the second term) reveals coupling between the 
charge(holon) and spin(spinon) degrees of freedom. Such 
coupling is ignored in other proposed theories [14-17]. 

Hubbard-Stratonovich transformations for the hop- 
ping, spinon pairing and holon pairing orders leads to 
the partition function, 



Z = j VfVbVxDA^ VA^ 



drC. 



(3) 



where Cef f = Co + Cf + Cb is the effective Lagrangian 
with 



J{l-xf 



<i,j> 



<i,3> 

for the order parameter Lagrangian and 

i,a 

_JJ1_^ E {xt,/l/.. + c.c.} 



(4) 



<i,j>,a 



J(l - x) 



E {^^UklfjA - fiAhl) + c-c.} (5) 



<t,3> 

for the spinon sector and 

Ci, = J2 biidr - li!')hi - t J2 {x:Mbj + cc.} 

~YH\'{^Vbib,+c.c.} (6) 

<i,j> 

for the holon sector. Here X; A-^ and A*" are the hopping, 

spinon pairing and holon pairing order parameters respec- 
tively, /i-' {/J') is the spinon (holon) chemical potential and x, 
the hole concentration. 

The SU(2) theory introduces the holon doublet 



hi = Qi 



b2i 



and the electron operator is then written. 



(7) 



(8) 



v2 

with the spinon doublet 'i'n = Qi ( ) = ( "(t^ ) and 
=gi( % ) = ( "^t't ) • Here gt is the SU(2) rota. 



/it / \ fil ^ 

tion matrix. Inserting Eq.(8) into the t-J Hamiltonian and 
following same procedure with U(l) case we obtain the effec- 
tive Lagrangian, 



jC^JP = Co + jCf + Lb, 



(9) 



where 



Co 



J(l - xf 



<ij> 



+ ^ E |Af,|^fe|A^,„,r+x^y (10) 

<i,j> \a,0 J 

Cs=Y,^\dr^i 
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<i,3> 
<i,3> 



Xij 



hi + C.C. 



(12) 



with >I>i = 



/.I 
ft 



It is quite encouraging to realize that 



the third term in Eqs.(6) and (12) manifests the composite 
nature of spinon pairing and holon pairing to allow for the 
formation of Cooper pairs. For the comparison of U(l) and 
SU(2) theory, we rewrite the holon hopping term(the second 
term in Eq.(12)) by using Eq.(7) as 



h 



Xij 



hi 



( bl bl ) 



Xij 

-A// 



-A: 



/ 

ij 
Xij 



bij 
b2j 



= b] 



.(1) 



-A^ 



/(I)* 



9 9j 



X'''ij bjbj, 



(13) 



where Xij"* and Af^'"'-' arc the U(l) mean field hopping and 
spinon pairing order parameters respectively. Thus the SU(2) 
order parameter matrix is related to the U(l) order parameter 



matrix U, 



(7(1) 



-A: 



/(I)* 



,ii) 



by gauge transformer 



tion, 



SU(2) 



Xij 

-A/; 



'A 



-Xij 



^ 9iUi^ 



c/(i) 



9]- 



(14) 



As a result, in the SU(2) theory the low energy order parame- 
ter fluctuations can be included and allows the appearance of 
holon doublet [16]. Such SU(2) treatment resulted in an arch 
shaped Tc line with the prediction of more realistic optimal 
doping rate than the U(l) theory [13]. 



B. OPTICAL CONDUCTIVITY IN U(l) AND 
SU(2) SLAVE-BOSON THEORIES 

We obtain the optical conductivity cr(a;) of an isotropic 2-D 
medium by evaluating the current response function Hxxi'jj), 



dJx{L0) 



1 d^F 



(15) 



Here Jx is the induced current in the x direction and Ex, the 
external electric field. F is the free energy, A^, the elec- 
tromagnetic field and w, the frequency of applied electric 
field. The total response function, Ilxa: = ^xx^ + ^xx'^ is 
the sum of the paramagnetic response function given by the 
current-current correlation function Hxx^ {v' — r,t' — t) =< 
jx(jc' ,t')jx{v,t) > - < jx{r',t') >< jx{r,t) > with the cur- 
rent operator jx{r,t) = ii(4, ^ ^(t)cr,cr(t) - £4,<r(i)cr+x,CT(t)) 



and the diamagnetic response function associated with the av- 
erage kinetic energy, Uxx^ =< K >= — t ^(4,<TCr+x,<T -|- 
H.C.) [20]. 

Recently, using the U(l) slave-boson theory we [18] studied 
the current response function contributed only from spinon 
singlet pair cxcitations(the amplitude fluctuations of the 
spinon pairing order parameter) and U(l) gauge fields(phase 
of the hopping order parameter) upto the second order, by 
using 



a,/3=a,|A/| 



n%„(n'' + n^)Jnv 



n^n'' 



+ ■ 



n 



+n-' 



aA nk+nf 



' W+TTT- 



(16) 



n 



where 11^(11'') is the spinon(holon) response function associ- 
ated with the internal and external gauge field (a and A); 

is the spinon(holon) re- 



n 



oA — dad\Al I y'^'^aA — dad\A} 

sponse function contributed from both the gauge fields and 
the spinon pairing field and H^^, H'aa and II^iAi the response 
function contributed from the spinon pairing field. The first 
term represents the loffe-Larkin rule [19] which results from 
the back-flow condition associated with gauge flcld fluctua- 
tions. On the other hand the second term is attributed to the 
spinon singlet pair excitations. It is noted that both terms in 
Eq.(16) contain the effects of the coupling between the charge 
and spin degrees of freedom. 

For completeness it is essential to include the effects of the 
phase fluctuations of hopping and spinon pairing order pa- 
rameters in addition to the amplitude fluctuations of hopping 
and spinon pairing order parameters. Thus we write the order 
parameter matrix. 



SU(2) 



-A^ii 

Ij^ 

"A''^ ij ~Xij 



X ij 



^/\xij? + WijV 



_ e cos 

-'nij\ ^g-i0ij si^Q.. _e'«« cos^ij 



e"-^*' sin 6*1 



(17) 



where r^ij = \/\Xij\^ + lA-*" ijP is the 'amplitude' of the order 
parameter matrix, aij{l3ij), the phase of the hopping (spinon 
pairing) order parameter and 9ij, the relative phase angle of 
the hopping and spinon pairing order parameters. It is noted 
that a, (3 and 9 represent the order parameter phase fluctu- 
ations in the sense that they are low energy fluctuations, r) 
represents the amplitude fluctuations of the order parameters. 

Allowing both the phase and amplitude fluctuations of the 
order parameters, we obtain the current response function. 



a.b=a.,!3 .0 .rj 



= tf- ^ n^„(n'' + n^):;nV, 



(18) 



a,b=0,ri 
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where = (u'^y = with X,Y = a,l3,e 

and rj. Here we used the fact that H^^^ = 11^^ = in the 
second hne. Therefore there is no coupUng between the exter- 
nal electromagnetic field A and the phase fluctuation modes 
a and /3. For details, we refer readers to Appendix A. 



III. COMPUTED RESULTS OF IN-PLANE 
OPTICAL CONDUCTIVITY 



3, 



e 

a 



3, 



0.1 
0.08 
0.06 
0.04 
0.02 




0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




1.5 

s 

a 



-S 0.5 



Fig. la; x=0.07 (Tc = 0.025t, T* = 0.062t) 
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Fig.lc; x=0.23 (Tc = 0.023, T* = 0.028) 
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FIG. 1. Computed optical conductivities as a function of 
temperature for x = 0.07(underdoped), x = 0.12(optimally 
doped) and x = 0.23(overdoped) cases with the antiferro- 
magnetic Heisenberg coupling strength of J=0.3t 

We present the predicted optical conductivities by allowing 
the order parameter fluctuations involving the parameters a, 
(3, 6 and rj upto second order in the framework of the SU(2) 
slave-boson theory of Lee and Salk [13] by choosing J=0.3t 



for the underdoped (a:: = 0.07), optimally dopod(a:: = 0.12) and 
overdoped(a; = 0.23) regions (Fig.l). Although not reported 
here, we find no qualitatively marked changes with the varia- 
tion of J in the range of O.lt < J < QAt. Here the bose con- 
densation temperature(To) and pseudogap temperature(T*) 
are given by our previous SU(2) results [13]. It is reminded 
that superconductivity is characterized by the order param- 
eters of holon pair and spinon singlet pair. The spinon sin- 
glet pairs appear at temperature below T*. Below the su- 
perconducting temperature Tc both the spinon singlet pair 
and holon pair order exist. The predicted results show the 
(Drude)Peak-dip-(MIR) hump structure at temperature be- 
low the pseudogap temperature T*(T < T*) and the hump 
disappears above T*{T* < T) indicating that the hump is 
originated from the spin singlet pair excitation, namely the 
antifcrromagnctic spin fluctuations of the shortest possible 
correlation lengths. Here we would like to point out that 
the predicted trend of both the shift of the hump position to 
lower frequency and the growth of the hump peak height with 
increasing hole concentration is consistent with observations 
[21,22]. However, quantitative results in hump position shift 
and hump peak height fail to agree with observations [3,4]. 



J=0.30t, x=0.07 (Tc = 0.025t, T* = 0.062t) 
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FIG. 2. Comparison of the optical conductivity 1) without 
fluctuations, with the fluctuations of 2) a, /3, 3) 9, 4) r] and 5) 
a, P, 0, rj; Result of case 1) and 2) is exactly same by virtue 
of n$i„ = n'Xp = 0(See Appendix A) 

To thoroughly examine the origin of the peak-dip-hump 
structure, below we discuss predictions made by systematic 
changes of various phase fluctuation modes of q , /3, and 9. 
Here a and (3 are the phases of hopping and spinon pairing 
order parameters respectively and 9, the relative phase angle 
between the two order parameters (see Eq.(17)). Fig. 2 above 
is a replot of Fig. la orfly for the case of T = O.OOlt. Without 
the inclusion of order parameter fluctuations only Drude peak 
appears. The same result is obtained with the allowance of 
only a and /3 fluctuations. Thus the computed optical con- 
ductivities with and without the phase fluctuations of a and 



P are identical. This is because fl^ 



n 



Af3 



■ owing to the 



SU(2) symmetry(see Appendix A for proof) as shown in Fig. 2 
above. The '9' mode fluctuations resulted in a broader Drude 
peak and some absorption at 200cm~^ < lj < 600cm~^. 
The 9 mode fluctuations represent an excitation mode as- 
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sociatcd with the amplitude fluctuations of both the hopping 
and spinon pairing order parameters which are massive exci- 
tations. The hump structure did not appear with the 9 mode 
fluctuations. Upto now we discussed the predicted optical 
conductivity associated with the phase fluctuations involved 
with a, f3 and 6. We now investigate the amplitude fluctua- 
tions o f the order p arameters involved with the '77' mode with 
r? = y^xP + |A7|2. shown in Fig.2 only the rj mode fluctu- 
ations give a peak-dip-hump structure in the optical conduc- 
tivity. However the width of the Drude peak is affected by the 
inclusion of the low energy excitations a, /3 and 6. The phase 
fluctuations cause further reduction in the width of Drude 
peak as compared to the case with the inclusion of only 77 fluc- 
tuations. Although not shown here, we observe that the phase 
fluctuations in the SU(2) cause a substantial reduction in the 
width of the Drude peak compared to the U(l) result but the 
position of the hump remains nearly identical. The SU(2) re- 
sult is in closer agreement with observations [3,4] compared 
to the U(l) case. On the other hand, the height of the hump 
is found to be lower in the SU(2) case as compared to the 
U(l) case. The hump structure is largely originated from the 
amplitude fluctuations of the 77 mode. Both the U(l) and 
SU(2) showed that the hump is originated from the antiferro- 
magnetic spin fluctuations of short range order, although we 
cannot deduce separate contributions from the two different 
amplitude(|x| and |A-^|) fluctuations in this SU(2) calcula- 
tion. This is because the amplitude fluctuations of both hop- 
ping order and spin (spinon) pairing order parameters occur 
in the background of antiferromagnetic fluctuations. 

For an additional check on the origin of the hump structure 
in a different aspect we computed the optical conductivity us- 
ing the Lanczos exact diagonalization method for a two hole 
doped 4x4 lattice for various values of Heisenberg antiferro- 
magnetic coupling strength J. The peak-dip-hump structure 
appears only for non-zero values of J and the position of the 
hump increases linearly with J. This strongly indicates that 
the hump is originated from the spin-spin correlations. This 
trend is in agreement with our present calculations. In Fig. 3 
we show the computed results of the hump position as a func- 
tion of the antiferromagnetic coupling strength J for compar- 
ison between the Lanczos exact diagonalization method (with 
one and two hole doped cases 4x4 square lattices) and U(l) 
and the SU(2) slave-boson theories. 
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FIG. 3. Antiferromagnetic coupling dependence of hump 
position for comparison with the exact diagonalization cal- 
culations with two hole for 4x4 lattice at T = QK and 
our slave-boson results at x = 0.07(SU(2)) and 0.05(U(1)) 
at T « Tc, that is T = O.OOlt 

We note that the peak locations of the hump obtained for 
both cases are sensitive to the variation of the antiferromag- 
netic coupling strength J, by showing a linear increase. This 
implies that the antiferromagnetic spin fluctuations are re- 
sponsible for the appearance of the hump structure. Thus we 
conclude that the hump structure in both cases is originated 
from the antiferromagnetic spin fluctuations, accompanying 
the amplitude fluctuations of the spinon pairing order pa- 
rameter at the hot spots. For comparison with observation 
we now present the predicted hump position with the varia- 
tion of hole concentration and temperature. In Fig. 4 the peak 
position of hump is seen to remain nearly constant as a func- 
tion of temperature far below T*; T* = 0.062f for x = 0.07, 
T* = Omt for X = 0.12 and T* = 0.028t for x = 0.23. 
In general, the predicted hump position tends to shift to a 
lower frequency with increasing hole concentration and with 
temperature, in agreement with the optical conductivity mea- 
surements [3,21,22]. As shown in Fig.l, at temperature close 
to the pseudogap temperature the hump structures are seen 
to too rapidly disappear with increasing temperature com- 
pared to observations. This implies that the antiferromag- 
netic spin fluctuations of short range order other than the 
spin singlet pair(excitations) are important even above T*. 
Such antiferromagnetic fluctuations of short range are con- 
sistent with NMR measurements [23] and other theoretical 
studies [7,10,11]. 
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IV. SUMMARY 

In the present study, we studied qualitative aspects of 
physics involved with the optical conductivity for the two- 
dimensional systems of strongly correlated electrons for a wide 
range of both hole doping(under-, optimal and over-doping) 
and temperature(T < Tc, Tc < T < T* , and T* < T) 
without resorting to any empirical parameters. By observing 
the linear increase of the hump position with the antiferro- 
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magnetic Hcisonborg coupling strength J, we showed that the 
hump structure in the optical conductivity is originated from 
the antiferromagnetic spin fluctuations which accompany the 
spin(spinon) singlet pair excitations(i.e., the antiferromag- 
netic fluctuations of the shortest possible correlation length) 
and that the AF fluctuations of larger correlation lengths arc 
also important for causing the hump structure. This obser- 
vation is consistent with NMR measurements [23] and other 
theoretical studies [7,10,11]. The shift of the hump position to 
lower frequency and the increase of the hump peak height with 
increasing hole concentration is sec to be consistent with ob- 
servations. We have demonstrated that the 9 mode phase fluc- 
tuations which represents the amplitude fluctuations of both 
the hopping and spinon pairing order parameters sharpens 
the Drude peak of the optical conductivity, yielding markedly 
better agreement with observations compared to the U(l) the- 
ory and that the hump structure is not affected by the a, (3 
and 6 mode fluctuations, that is, the phase fluctuations in- 
volved with the hopping and spin singlet pairing order pa- 
rameters. As pointed out earlier, our computed results fail 
to agree with observations in the magnitude of hump height. 
To achieve quantitative accuracy in the future study it is de- 
sirable to include the next nearest hopping term(t' term) in 
the t-J Hamiltonian and the effects of the AF fluctuations of 
larger correlation lengths. 
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APPENDIX A: PROOF OF = U% = IN THE 
SU(2) THEORY 

From Eq.(12), the holon action in the presence of the ex- 
ternal electromagnetic field is written. 



rl3 

Sb{h,a,P,e)= / dr^/ij(ar 
Jo i 

<i,j> ^ ' 

-{^^ E ( J -1 ) K')" + ' (Ai) 



<i,j> 



where we used Xij = Vij^^"'^ cosOi and ij = rjijC^^'i sinOij . 
Under the local SU(2) transformation holon field and the order 
parameters are transformed as 



U'ij = giUijS^., 



(A2) 
(A3) 



where gi = e*2 '®' is the SU(2) transformation matrix 



and Ui; 



'■^ 1 the order parameter ma- 



ij Xij 



trix. Here f are Pauli spin matrices and Qi, the rota- 
tion angle. Considering the SU(2) rotation of ±7r around 
the 2-axis the SU(2) transformation matrix is given by 



1 



and Uii is transformed as 



e 'JcosS;. -e 'JsinS., 



\ — e »3 sine', — e '3 cosS... 

-( ^\ix+iv ( l\ ( cos $1 -e'^isinSiN 

- I V -1 ; -e-'"! sine, -e*°i cos 6, J 



*^cos^i — e sin 01 
e''^'- sin 5; — e~*"' cos^; 



(A4) 



where J = i-|-a;orj =i + y and (— 1)*^+*« x (— = 
-1. Thus we obtain = -a^-, P'^j = -Pij, O'^j = 9ij 
and r]lj = rnj . Then the holon action is given by 



St(h,A,a,0,e,ri) 



He'-*..„..fttJ„, /e-'°'^ cose., -e'^-^ sine.A t„.ft . 

<i,j> 

<i,j> 

I dr'^h'lid.r - n)h'i 

■ 

^ ^-^ y V-e 'JsinS'.-e 'J cos 61'./ J 



<ij> 



: S''{h',A,-a,-/3,e,ri). 



(A5) 



Thus we find that the holon free energy is the even func- 
tion of a and /3: 



--In / m'e^^''^''''^'-"'-^'^''') 



= F\A,-a,-(3,e,7^), 



(A6) 



where we used the fact that the holon free energy is in- 
variant under the local SU(2) transformation. It is then 

obvious that the holon current response function gj^^ 
satisfies 



Aa 



dAda 



a=l3=Q 



d^F\A,-a,-l3,e,ii) 



dAda 



a=/3=0 



(A7) 
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Thus we obtain = 0. Similarly = 0. The 

physical meaning of 11^^ = 11^^ = is that there is 
no coupling between the external electromagnetic field 
A and the order parameter phase fluctuations a and (3. 
Thus a and (3 mode fluctuations are not affected by the 
external electromagnetic field and the external electro- 
magnetic field can not screen the fluctuations of these 
modes. 
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